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HIGHLIGHTS 


• A wood pellet steam engine micro CHP device has been tested. 

• A micro CHP data driven dynamic model has been implemented in TRNSYS environment. 

• The model takes into account transient and steady state phases. 

• The model is suitable with building energy simulation. 
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A wood pellet micro combined heat and power device (pCHP) has been tested in order to characterize its 
performances in steady and transient states. A dynamic model based on these experimental in¬ 
vestigations has been developed in order to predict its energy performances and its pollutant emissions. 
The model is designed with a few parameters experimentally accessible. This model has been imple¬ 
mented in TRNSYS numerical environment. This work focuses on the experimental investigations and on 
the model description. The modelling approach is based on a physical part (an energy balance on the 
entire device and a combustion model), and on an empirical part (correlations for the fuel power input 
and for the thermal and electrical outputs). The model characterizes the pCHP behaviour for different 
part load ratios (PLR) (power modulation). The dynamic phases with start-up and cooling phases are also 
taken into account. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Micro combined heat and power (pCHP) is an emerging means 
in order to produce simultaneously heat and power with high ef¬ 
ficiency. CHP (or cogeneration) valorises the residual heat gener¬ 
ated during power production by means of thermodynamic cycles 
such as Rankine or Stirling cycles. This heat production is used for 
space heating and domestic hot water (DHW) while the electricity 
production can be used in building or exported on the grid 
compared to centralized power plants such as nuclear, gas or coal 
power plants that use cooling towers to remove the heat. CHP 
power plants increase the energetic efficiency and decrease the 
environmental impacts by consuming less primary energy 1-3 . 
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CHP is already widespread in industrial applications or building 
applications in district heating for a power range of MW in order to 
be suitable for electrical and heating needs. pCHP corresponds to 
power compatible with residential or commercial heating and 
power needs. pCHP systems have a power output below 20 l<W e i 
[4]. The technologies which are used are different from centralized 
power plants [5] and the electrical efficiencies are lower. These 
decentralized systems also reduce losses on distribution networks 
[6]. The main conversion processes which are used for pCHP ap¬ 
plications are: internal combustion engines (ICE), micro turbines, 
Stirling engines, Ericsson engines, steam engines (SE) and fuel cells 
[1 , [7 . Steam engines present certain advantages compared to the 
other technologies. The systems are simple by using free pistons or 
rotary engines, their design is cheap and the steam is easy to 
generate with any fuels such as gas, oil, solar energy or biomass [8]. 
At micro scale, the use of the Rankine cycle thresholds power ef¬ 
ficiency at about 10% but heat recovery can be very high [8]. Dong 
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et al. [9] made a review of all technologies suitable with biomass 
pCHP applications. They conclude that the biomass fuel can be 
applied to all classical pCHP technologies listed before and they 
describe all the energy conversion processes of the biomass: 
combustion for steam generation, gasification, pyrolysis or bio¬ 
chemical processes. 

A power modulating wood pellet steam engine pCHP (WP SE 
pCHP) device has been tested to characterize its behaviour in steady 
and transient states. The objective of this paper is to describe the 
experimental characterization of a WP SE pCHP leading to a per¬ 
formance based model dedicated to annual building simulations. 

The behaviour and the performances of pCHP systems are still 
little known; in particular WP pCHP devices are little studied 10]. 
Among the existing commercial pCHP devices, only some of them 
use renewable fuels. Table 1 sums up the existing devices on the 
market or under development. 

Some works have been achieved on biomass pCHP devices. 
Thiers et al. [10] propose a study on a commercial WP pCHP device 
using a Stirling engine (Sunmachine). Thermal and electrical out¬ 
puts have been studied. Some works propose to couple a classical 
boiler (using gas, fuel or biomass) with an organic Rankine cycle 
(ORC). Cardozo et al. 11 connected an overfed WP burner to a 
gamma type Stirling engine. This configuration led to an overall 
efficiency of 77% and an electrical efficiency of 9%. Guoquan et al. 
[12] underline the relevance to connect a biomass boiler used as 
steam generator to a steam engine. They carried out a prototype 
which links an ORC engine with a biomass boiler but the electrical 
efficiency only reaches 1.48%. Cordiner et al. 13] used a compressor 


Table 1 

Review of commercial WP (or which accepts any fuel) pCHP devices. 


Manufacturer 

Apparatus 

Power 


LHV efficiency 



Electric 

[kW el ] 

Thermal 

[kW th ] 

Electric 

[%] 

Thermal 

[%] 

Stirling engine 

0 OkoFEN 

Pellematic 

Smart_e 

1.1 

15.8 

7 

95 

Sunmachine 

Pellet 

1.5-3 

4.5-10.5 

20-25 

65-70 

Rankine: scroll turbine (module to connect to external heating source) 


Enerftech 

EnefcogenPLUS 

5 

33 

12.5 

82.5 

Eneftech 

Enefcogen 

GREEN 

5-10 

45-80 

10.5 

84.2 

Rankine: scroll turbine 





Novotek 

industry 

COGEMAX 14 

0.6-1.5 

4.8-12 

10 

82 

Rankine: steam engine (linear free piston) 




Exoes/Otag 

BISON 

0.3-1.6 

3-16 

9 

83 

Rankine: steam engine (module to connect to external heating source) 


Exoes 

SHAPE 

3 

50 

6 

84 


ready. 

introducing phase. 

testing and optimization phase. 


as steam expander connected to a biomass furnace. Experimental 
studies led to an electrical efficiency of 13% with an electrical 
output of 70 kW. Alanne et al. [14] studied the implementation of a 
thermoelectric cogeneration system inside a biomass boiler. 
Theoretically, electrical efficiencies of 9% can be reached. They also 
studied the coupling between a rotary steam engine and a biomass 
boiler where electrical efficiencies can also theoretically reach 9% 
[8]. Creyx et al. [15] studied the opportunity to couple a hot air 
engine based on an Ericsson cycle with a biomass boiler where 
theoretical thermodynamic efficiencies of about 37.6% can be 
reached. Prando and al [16] propose to couple a biomass gasifica¬ 
tion reactor to an internal combustion engine for residential ap¬ 
plications. Finally, Angrisani and al [17] developed a concept of 
solar-biomass pCHP system based on a fluidized bed and a Stir¬ 
ling engine designed to produce about 1 kW e i. 

Different modelling approaches exist to simulate the behaviour 
of a pCHP device. A full physical model requires the modelling of 
the combustion reaction, the thermodynamic cycle and mechanical 
and electrical conversion with non-ideal efficiencies. This approach 
needs data on the engine and the model supposes low computing 
time steps which are not suitable with annual dynamic building 
energy simulations. The usual time steps for accurate thermody¬ 
namic models (from the order of 10 3 to 10~ 6 s) are very different in 
comparison to buildings (from 1 min to 1 h). Moreover, to correctly 
model an engine, many features have to be known like heat 
exchanger surfaces, work pressures, geometry, etc, which are 
generally unknown. Alanne et al. [8 implemented a thermody¬ 
namic model of a rotary steam engine using a Rankine cycle by 
taking into account work pressures, isentropic efficiency or 
expansion volumes. Creyx et al. 15] developed a thermodynamic 
model on expansion and compression phases of an Ericsson engine 
coupled to a biomass boiler for design optimization applications on 
the operating pressures and temperatures. This model is not used 
for energy building simulations. Lontsi et al. [18] developed a 
thermodynamic model of a Joule cycle reciprocating Ericsson en¬ 
gine for pCHP systems which takes into account start-up phases. 
This system would be suitable with renewable fuels such as 
biomass fuel. The simulation time step is 2.10 -3 s, and the model is 
used to optimize the design of the engine. 

Semi-physical simplified models of internal combustion engine 
and Stirling engine have been developed for annual energy building 
simulation. In the International Energy Agency (IEA)/Energy Con¬ 
servation in Buildings and Community Systems (ECBCS) Annex 42, 
Beausoleil-Morrison et al. [19] propose a general physical and 
parametric model without thermodynamic considerations which is 
suitable with biomass applications. The Annex 42 gives a literature 
review on models of mini-and pCHP tools for dynamic simulation of 
buildings. The Annex 42 proposes a method to model pCHP dis¬ 
tinguishing steady states and transient phases through grey box 
models. It distinguishes 4 phases: stand-by, warm-up, steady state 
and cool down phases. These models use thermal inertia to model 
transient phases. Steady state performances are calculated by using 
correlations based on the main parameters that influence thermal 
and power inputs and outputs: the cooling water temperature and 
mass flow. The required time step is suitable for annual dynamic 
thermal simulation (a few minutes). Lombardi et al. [20] have 
simplified this model by implementing a linear difference model 
based on linear coefficients and reference parameters which are 
determined by a few experimental tests. 

Furthermore, “black box” models have been developed. Thiers 
et al. 10] propose a simple model on a WP SE pCHP by using a fully 
empirical modelling approaches based on experimental in¬ 
vestigations. They have fitted a quasi-static model from a mea¬ 
surement campaign conducted on the Sunmachine. The pre¬ 
heating and cooling phases are characterized only on the 
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combustion chamber temperature. This model is strictly limited to 
the tested product and to the achieved test conditions (source 
temperatures, water flow rates, fuel characteristics, heating loads). 
The thermal and electrical outputs are calculated only in steady 
states according to the inlet cooling water temperature. 

Indeed, Conroy et al. 21 ] made a data driven numerical model of 
a Stirling engine pCHP device based on field tests. The accuracy of 
the measurements is not given and the fitting only depends on time 
without physical meaning. 

Our objective is to build a simple and accurate dynamic model 
which requires the fewest parameters. These parameters have to be 
easily reachable by experimental investigations or manufacturer 
data in order to carry out optimization simulations on annual en¬ 
ergy building simulations at low simulation time steps (up to 
1 min). The parametric model developed in the Annex 42 is the 
closest to these objectives. 

Many studies are carried out in order to assess the primary 
energy gains or the environment impact of pCHP without taking 
into account the transient behaviour of the pCHP devices and by 
only using constant heat to power ratio or constant efficiencies 
[1,6,22,23 . Hawkes et al. [24] investigated the influence of the 
simulation time step of the electrical and thermal demands on the 
energetic, economic and environmental performances of pCHP in¬ 
tegrated to buildings. They conclude that the results are very sen¬ 
sitive to the time step when varying between 1 h and 5 min: low 
temporal resolution lead to pCHP oversizing, economic gains un¬ 
derestimations and environmental impact overestimations. Only a 
power ramp rate, the minimum/maximum power output and the 
efficiencies in steady state are used for the pCHP model. The tran¬ 
sient behaviour of the pCHP devices is not taken into account. The 
objective of the developed model is to consider the dynamic 
behaviour of the pCHP in order to assess the performances of pCHP 
integrated in residential buildings. 

This paper focuses on the test bench description, on the 
experimental investigations and results of a WP pCHP using a linear 
free piston steam engine and on the model description. 

2. Experimental study 

2.1. Description of the wood pellet gCHP test bench 

The experiments have been carried out on the WP pCHP BISON 
from EXOES. The WP pCHP unit has a steam engine (linear free 
piston) supplied by a steam generator which can produce steam up 
to 20 bars and 400 °C. According to the manufacturer, the thermal 
output can vary from 3 to 16 l<W t h and the electrical output can vary 
from 0.3 to 1.6 kW e i. The device controls the incoming combustion 
air flow to vary the load. The supply of pellets is controlled to 
maintain a constant level of pellet inside the combustion chamber. 
The pellet supplier is indirectly controlled by the device to control 
the part load ratio (PLR). The test bench has been designed in order 


to study instantaneous heating or full heat storage configurations 
as is shown in Fig. 1. For the heat storage configuration, the PLR is 
controlled according to the tank top temperature. For the experi¬ 
mental tests, the volume flow rate and the PLR are fixed and the 
inlet water temperature is controlled according to the set point by 
playing on a 2 ways valve on the heating circuit. 

The electrical power output and electrical auxiliary consump¬ 
tion (starter, pumps, fan, WP feeder, air compressor, and ash 
extractor) are measured by two separate watt meters. Thermal 
output is assessed by using a flow meter and two temperature 
probes. A gas analyser is used to perform measurements on the 
exhaust gas. The fuel mass flow is assessed by using a scale on 
which the storage tank is placed. The WP feeder is mounted on a 
support in order to avoid pressing on the WP storage tank. The 
connection between the feeder and the WP storage tank has been 
designed in order to avoid any contact. 

The pCHP unit is composed of 4 levels (Fig. 2). The first one at the 
top is composed of the steam generator where combustion occurs. 
There is also an air compressor and a high pressure bottle in order to 
clean up the pipes of the heat exchanger (steam generator) each 
1000 s. The second one in the middle is composed of the steam 
engine which is supplied by steam coming from thin highly isolated 
pipes. Under the combustion chamber, there is a starter which heats 
up the WP inside the combustion chamber to launch the combus¬ 
tion. The third one is composed of the tank which receives the hot 
condensates from the steam engine. There is also a high pressure 
pump (HP pump) which is used during the start-up and cooling 
down phases. It pumps the low pressure water from the heat 
exchanger to supply the steam generator at high pressure. Hot 
condensates from overheated steam are stored in the tank before 
being pumped up to the heat exchanger which is on the last level at 
the bottom. There are also 3 pumps: one for the primary heat circuit, 
one for the heat distribution circuit in the building and one which is a 
vacuum pump to control the low pressure level in the tank. 

2.2. Metrology 

Table 2 summarizes all the metrological means implemented on 
the test bench. The global uncertainties for each quantity are 
calculated according to the law of propagation of uncertainty 

[25-28]. 

2.2.2. Fuel power 

The fuel power which is supplied to the device is calculated by 
Eq. (1) based on the fuel mass flow and on the LHV which is derived 
from laboratory tests on WP samples. 

Pfuel = m fuel LHV (1) 

The fuel mass flow is calculated with a scale which has been 
chosen in order to minimize the uncertainty. Since the greater the 


Electricity 

consumed 

Electricity 

generated 

Wood 

pellets 


1 
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W 3 way valve 
( ^T >/ Safcty valve 
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^Expansion tank 
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Y Backflow preventer 

M Non return valve 
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Fig. 1 . Operating scheme of the test bench (with the control strategy principle on the tank and inlet temperatures). 


































1042 


J.-B. Bouvenot et al. / Applied Thermal Engineering 73 (2014) 1039—1052 



Fig. 2. Technical description of the BISON (1: WP supplier, 2: exhaust gas, 3: fan, 4: steam generator, 5: ash extractor, 6: linear generator, 7: high pressure pump, 8: tank, 9: filter, 10: 
heat exchanger, 11: primary circuit pump, 12: buffer tank, 13: compressed air tank, 14: starter, 15: secondary circuit pump). 


Table 2 

Metrological means and associated uncertainties. 


Type 

Parameter 

Notation 

Metrologic means 

Uncertainty 

Fuel input 

Consumed fuel 

fftfuel 

Scale (Mettler Toledo ICS425) 

±10 g 


LHV 

LHV 

DIN EN 14,774-2 

±0.1 kWh/kg 


Fuel and air temperature 

TmbTfuel 

4 wires Ptl 00 

±0.2 K + 0.05% 


Fuel power 

Pfuel 

Calculated Eq. (2) 

±3% 

Heating output 

Water temperature 

Tw,i> Tzw.o 

4 wires Ptl 00 

±0.2 K + 0.05% 


Volume flow 

l^cw 

Blade flow meter 
(Biirkert Ml2) 

±0.15 lL/mn + 2.5% 


Thermal power 

Qhx 

Calculated Eq. (1) 

±4.5 to 7.3 % 
(according to PLR) 

Electric output 

Gross electric power 

P gross 

Energymeter 
(Socomec DIRIS A20) 

1% 

Electric input 

Consumed power 

Paux 

Energymeter 
(Socomec DIRIS A20) 

1% 

Tank temperature 

Storage tank temperature 

Tank (x4) 

4 wires Ptl 00 

±0.2 K + 0.05% 

Exhaust gas 

Temperature 

0 2 , NO x , CO, C0 2 


Analyser ECOM J2KN 

— 


weight range is, the greater is the measurement uncertainty, a 
compromise has been found in order to threshold the WP storage 
tank capacity to a 8 h full load test. The mass flow is smoothed on 
the whole test so as to avoid noise on the measurements linked to 
the feeder vibrations in the WP mass and WP dynamic movements 
inside the tank (Fig. 3). The mass flow corresponds to the slope of 
the mass variation. The uncertainty is mastered by choosing the 
test time according to the PLR in order to reach an uncertainty 
below 1% (Table 2). 


Pth = 


Qhx 

Auei 



2.2.3. Electrical power 

The electrical power is recorded by 2 W m, one measuring the 
incoming electricity for the operation of the auxiliary devices and 
the other measuring the gross electricity which is produced by the 


2.2.2. Thermal power 

Thermal power is calculated by using the current equation of 
heat transport (Eq. (2)). 

Q.HX — p(^cw,i) '^cW (O:w(7cw,o) '7cw,o — ^cw(^cw, 0-W (2) 

The variations of the thermal properties with the temperature of 
the heating water are taken from Holman 29 . The density is 
calculated with the cooling water inlet temperature because the 
flow meter is located on the backflow, just before the temperature 
probe. The thermal efficiency is given by Eq. (3): 



time [s] 
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steam engine. The uncertainty is below 1% ( fable 2). The gross 
electrical efficiency is given by Eq. (4): 


Vel = 


f^gross 

^fuel 



We also define the power to heat ratio and the global efficiency 
by Eqs. (5) and (6): 


a = 


^gross 

Qhx 



V g = Vth + Vel 



3. Experimental results 

The experimental investigations aim at assessing accurate en¬ 
ergy performances of the device and well characterize the transient 
phases: the start-up and the cooling phases. 

3A. Steady state phases 

Experimental tests have been made on a PLR range of [22 %— 
98%] by step of about 10%, on a water inlet temperature range of 
[30 °C-45 °C] by step of 5 °C and on a water volume flow rate range 
of [5 L/min-15 L/min] with a step of 5 L/min. Tests made from 40% 
to 80% of PLR show the low influence of the water inlet tempera¬ 
tures and the flow rates on the thermal, electrical and fuel power 
(Figs. 4 and 5). The influence of these parameters is lower than the 
uncertainty or repeatability ranges of the measurements. The 
thermal efficiency is not impacted by the variation of the water 
flow rate of the heating circuit. Heat recovery and steam generation 
are separated. 

The variation range of the cold source temperature is very low 
(maximum gap of 15 °C) compared to the difference between the 
hot source (inside the combustion chamber up to about 900 °C) and 
the cold source (condensate steam in the tank) that makes an 
insignificant influence on the thermal and electrical performances. 
This device does not take profit of the exhaust gas condensation so 
the performances are low dependent on these 2 parameters. For 
any cooling water flow rate or temperature, the steam is produced 
by playing on the incoming combustion air and WP feeder. This 
steam is injected into the engine to produce power. Then, the steam 
is condensed in the buffer tank. The temperature of the buffer tank 





Fig. 4. Steady state results on fuel, thermal and electrical powers. 



20 30 40 50 60 70 80 90 100 PLR [%] 


Fig. 5. Steady state results on the global, thermal and electrical efficiencies and on the 
power to heat ratio. 


depends on the withdrawal of heat toward the heating circuit. 
There is no correlation between the temperature level in the tank 
and the physical properties (pressure and temperature) of the 
steam injected into the engine (which is controlled). During the 
tests, the volume flow and inlet cooling water temperature have 
been chosen to minimize the uncertainty range. The thermal, fuel, 
electrical and auxiliary powers are linear according to the PLR 
which is an internal parameter of the device (Fig. 4). These 4 
quantities do not tend toward 0 when the PLR tends to 0, so the 
validity range of the PLR for the linear correlations is only between 
22 and 98%. 

The efficiencies are highest for full load operations and the 
global efficiency can reach 95% on LHV (Fig. 5). The efficiencies 
decrease slightly with the PLR reduction up to 22%. The power to 
heat ratio is relatively stable around 0.1 up to 41% of PLR. This de¬ 
vice has similar energy performances to others existing or devel¬ 
oping biomass Rankine pCHP devices ( Table 1). 

Experimental investigations also show other main results that 
are useful for the numerical model such as the oxygen rate in the 
wet flue gas which is controlled at a constant value of 6.2% or the 
ash temperature at the bottom of the combustion chamber which is 
constant for any PLR at about 230 °C. 

The exhaust gas has also been studied. The CO and NO contents 
have been measured using a gas analyser (J2KN from ECOM) as well 
as the temperature during 2 h long tests. Fig. 6 gives the mean 
values of the CO and the NO volumetric contents according to the 
PLR. The NO production increases with the PLR although CO 
emissions are almost constant according to the PLR ranging be¬ 
tween 20 and 40 ppm. These values are very lower than values 
fixed by the European standard on wood pellet boilers (EN 12809). 
The NO emissions increase with the exhaust gas temperature which 
is linked to the PLR. The exhaust gas temperature will let to 
calculate the exhaust gas heat losses. As it has been noticed in Fig. 4, 
the exhaust gas temperature and the CO and NO volumetric con¬ 
tents do not tend toward 0 when the PLR tends to 0. The linear 
correlations for these quantities are again valid only for a PLR range 
between 22 and 98% which are the minimum and maximum PLR of 
the tested device. 


3.2. Start-up phase 

The pCHP device at the beginning uses the 2 pumps of the 
secondary and primary circuits, the fan and the electronic card. 
Then, after 20 s on average, the starter begins to warm up the 
combustion chamber with a mean electrical power of 1250 W for 
an average time of 300 s. The thermal production begins on average 
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Fig. 6. Steady state CO and NO contents in exhaust gas. 


after 450 s, the engine is bypassed and the steam supplies directly a 
coaxial heat exchanger on the heating circuit. The PLR set point is 
constant at 27% for each start-up phase. Then, the high pressure 
pump is triggered until the steam engine starts after 670 s on 
average. Fig. 7 gives the behaviour of each quantity during the start¬ 
up phase. The engine is used in steady state operation to pump 
water from the low pressure to the high pressure circuits by using a 
part of the mechanical energy of the engine. On average, the start¬ 
up cycle involves additional 142 Wh e i of electrical energy which is 
consumed by the start-up devices (starter, over use of the WP 
feeder and the high pressure pump) before launching the electrical 
production. In Fig. 7, the abrupt loss of thermal power after 700 s 
comes from the switch in the hydraulic circuit and the beginning of 
the conversion of a part of the heat in power. The total duration 
start-up cycle is on average of 1100 s to reach the steady state 
values of each quantity. Indeed, as soon as the starting cycle of the 
engine is triggered, the overheated steam is injected in the main 
cold circuit, that is why the system needs time to extract the heat 
from the system because the condensates tank need to be ther¬ 
mally loaded. The trend curves in Fig. 7 will be used to model the 
start-up phases. 

The fuel mass flow rate variation is also measured. There is an 
ignition in order to supply the combustion chamber by about 100 g 
of biomass fuel. Then, as soon as the starter turns off, the combustion 
starts and the feeder works for about 150 s in order to completely fill 
the combustion chamber with about 500 g of biomass fuel (Fig. 7 
between 350 and 500 s). Then, the mass variation is linear during 
steady state (Fig. 3) and the variation only depends on the PLR. The 
WP feeder operates according to the PLR to fulfil the combustion 
chamber in order to maintain the WP level as constant thanks to a 
light sensor. The higher the PLR is, the greater is the electrical con¬ 
sumption of the feeder. This increasing is integrated in the whole 
electrical consumption of the auxiliary devices (Fig. 4). 

3.3. Cooling phase 



Fig. 8. Cooling phase: thermal and electrical production and auxiliary consumption for 
PLR = 80%. 

occurring until the end of the combustion of the remaining biomass 
fuel. As soon as the engine is switched off, the high pressure pump 
is turned on intermittently during the whole of the cooling phase. 
The cooling phase ends when the high pressure pump stops, after 
1800 s on average. The auxiliary electrical overconsumption linked 
to this pump reaches on average 38 Wh e i. The trend curves in Fig. 8 
will be used to model the cooling phase. 

3.4. Modulating phase 

When the PLR varies, the device has a response time. There is a 
power ramp rate which cannot be exceeded to reach the PLR set 
point. Tests have been carried out to assess this speed which rea¬ 
ches 6 Wth/s for the thermal power and 0.6 W e i/s for the electrical 
production (Fig. 9). 

3.5. Exhaust gas 

The start-up and the cooling phases have been studied in order 
to assess the pollutant emissions which are higher during these 2 
phases. NO emissions are stable during the start-up phase although 
CO emissions are higher than steady state emissions (Fig. 10). 
During the cooling phase, both the CO and NO emissions are higher 
than steady state emissions. During the start-up phase, the CO 
emissions are higher because at the beginning, the combustion 
chamber is cold and the low temperatures slow down the oxidation 
of the biomass carbon and generates more CO [30]. During the 
cooling phase, an oversupply of air occurs that also cools down the 
combustion chamber. The NO emissions remain stable for each 
phase. 


During the cooling phase, any fuel is consumed and the thermal 
and electrical power decrease (Fig. 8). Electrical production de¬ 
creases up to a threshold of 500 W. The thermal production keeps 


10,000 

a> 

■t 

© 8,000 

a 

«3 £ 6,000 

f “ 

1 4,000 

a> 

£ 2,000 


starter 


R 2 = 0,9405 


QHX 


2,500 


-QHX trend 

’ pump + WP feeder .. 

( R 2 = 0,9405 - Pgross 

- Pgross trend 

- Paux 




0 100 200 300 400 500 600 700 800 900 1,000 1,100 time [s] 


Fig. 7. Start-up phase: electrical consumption, electrical and heating productions for 
PLR = 27%. 


3.6. Ash production 

The ash content has been assessed by recovering the ash during 
long tests (at least 24 h) in order to increase the accuracy. An ash 
content in the WP of 0.3% has been found on average. 



Fig. 9. PLR variation speed in response of a PLR step change from 27 to 90%. 
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Fig. 10. CO and NO emissions during the transient phases. 


4. Wood pellet micro CHP numerical model 

The developed model is mainly based on the previous experi¬ 
mental studies. The 4 phases described in the Annex 42 are studied 
separately. 3 parameters have been studied: the PLR, the cooling 
water mass flow and the inlet water temperature. The experimental 
results show that only the PLR plays on the performances. The main 
reason is that the studied device does not recover the latent heat 
from the condensation of the steam inside the exhaust gas. 

A model has been developed in the TRNSYS environment under 
the name of Type 254 in order to well characterize the physical 
behaviour of the pCHP. The model tries to match the objectives of 
simplicity and accuracy to be suitable with annual dynamic thermal 
building simulations this is the reason why it is designed as a grey 
box model with a physical part and an empirical part. The energy 
balance and the combustion model are physical and the main in¬ 
puts and outputs are calculated by using correlations according to 
the PLR coefficient. In steady state, the energy balance has to be 
respected: 


FFfuel + ^air + Pfuel — Fgross + Oloss + Q.HX + ^exh (7) 

In steady state conditions, the electricity power, the fuel mass 
flow and the heating power are assessed by correlations dependent 
on the PLR by using the linear differences model of Lombardi [20]. 
The combustion model is based on classical chemical reactions 
involving the combustion air, the fuel and the reaction products. 


CW ,0 


Qhx 

Ccw^cw 


+ T, 


CW,l 


(13) 


The Appendixes 1 and 2 present the transient states transition 
conditions and the algorithm chart with the main inputs and out¬ 
puts with internal calculations. 


4.2. Transient phases 

In addition to the steady state characterization, the start-up and 
cooling transient phases need to be modelled. These phases are 
modelled by exponential laws and delay times (Figs. 7 and 8 ). The 
modification of the PLR set point is not instantaneous in reality. In 
the numerical model a power ramp rate is implemented (Fig. 9). If a 
PLR set point occurs, the process value of PLR will vary according to 
the power ramp rate. For the electrical consumption of the auxil¬ 
iaries, the overconsumption of the electrical energy linked to these 
devices for each phase is divided by the transient phase duration: 
670 s for the start-up phase and 1800 s for the cooling phase (Eq. 
(14) and (15)). 

if t < fstart + Afstart = Faux (PLRstart) + . ^p Yt 04) 

^‘■start 

E P 

if t> t s top + Atf top P aux P = Paux (PLRstop) H q —p 

A f stop — A f stop 

(15) 

The Appendix 3 presents the modelling representation of the 
start-up and the cooling phases. 


4.3. Combustion modelling 

From the fuel mass flow rate, the fuel composition and the 
combustion air excess, typical combustion calculations are under¬ 
taken in order to assess the CO 2 mass flow rate and exhaust gas heat 
losses. 


4.1. Steady state phase 

The expressions of the delivered thermal output, fuel input, and 
gross electric power (in W), are defined as a function of the nominal 
conditions for each value with PLR n0 m = 100% (PLR <e [22 - 100]) by 
Eqs. (8)—(10) (LHV is given in kWh/kgf ue i) 

P f U ei = 1000- (a(PLR - PLR nom ) + n'CT jLHV (8) 

• nom 

Qhx = b( PLR - PLR nom ) + Q HX (9) 

Pgross — C(PLR - PLRnom) + PgSss (10) 

The net electrical power corresponds to the usable part of 
electricity produced by the engine. The other part of the production 
is self-consumed by the auxiliaries (pumps, fans, electronic card, 
etc.). The net electrical power is defined by Eq. (11). 


4.3A. Chemical composition and energy density of the wood pellet 

The model requires the mass fraction of each of the following 
elements on dry fuel: oxygen (O), carbon (C), hydrogen (H), the ash 
content and the water content. Either the chemical composition is 
derived from laboratory tests as has been carried out here ( Table 3), 
or default values are proposed: the chemical composition of the 
pellet is quite stable and the default values can be used with a good 
accuracy [31,32 . For the ash content, laboratory tests can be carried 
out or the standards [32] or the manufacturer data can be used. 

The LHV and the water content in the fuel are two key values of 
the model. Either they are also derived from laboratory tests or 
default values are proposed. Bernard [30] proposes a simple for¬ 
mula based on the mass composition of the dry WP and which gives 
the LHV value of the dry part of the WP (Eq. (16)). The LHV on the 
dry part takes into account the required energy for the trans¬ 
formation of the hydrogen in the WP and which is turned into 
water. 



gross 


aux 


(ii) 


Paux = d(PLR - PLR nom ) + P£T (12) 

For the model identification, the LHV of the fuel should be 
known or measured. The outlet cooling water temperature is 
calculated using the heat transport equation: 


Table 3 

Experimental and theoretical WP characterization. 



w 

^ash 

Y C dry 

Y H dry 

Y 0 dry 

LHV dry 

LHV 


[%] 

[%] 

[%] 

[%] 

[%] 

[kWh kg -1 ] 

[kWh kg 1 ] 

Experimental 

6.6 

0.34 

51.3 

6.2 

42 

5.29 

4.9 

Theoretical 

8 

0.4 

50 

6 

44 

5.34 

4.9 
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LHV = (1 - w) (34,030 Y c dry + 121,640 Y H dry 

- 12,540 Y odry ) -wL v 


(16) 


Here, the LHV, the chemical composition and the water content 
have been tested in a laboratory by using specific test standards 

(Table 2). 


4.3.2. Wood pellet combustion reactions model 

The chemical reaction of the dry part of the input fuel mass flow 
rate, without the ash part which is neglected, gives: 

C*HyO z(s) + (* + T - 7 ) (l + e air) 02(g) -^*C0 2 (g) 

y / y z\ 07) 

+ J H 2°(g) + ( x + 4 _ 2j e air°2 (g) 

Different values are calculated in order to know each mass flow 
rate of each component. 


rfifuei 


^fuel 

LHV 


(18) 


niair = nifuel(l -W) 


y _i_ y tl. 

A + 4 2 

M C X H y O Z 


(1 T" ^air) 


( Xo 2 Mo 2 + XN 2 Mn 2 ) 


XO , 


(19) 


The content of water in exhaust gas is composed from 2 con¬ 
tributions: the part which comes from the hydrogen content of the 
dry WP and the part which comes from the water content of the 
wet WP. 


The CO and the NO masses have been integrated during the 
transient phases in order to be used in the numerical model. For 
each starting cycle, an emission value for the CO and NO emissions 
is incremented like Persson et al. did for WP boilers and stoves [33]. 
Table 4 summarizes these values. During the steady state phase, the 
model takes into account the CO and the NO emissions by Eqs. (28) 
and (30). The gas analyser gives the CO and the NO volumetric 
content in the wet exhaust gas. Eq. (27) gives the volume flow of the 
exhaust gas and Eq. (29) gives the NO emission variation according 
to the PLR (Fig. 6). 


V ex h = ^fuel^i 


m 


W (1 - W) 


Mh 2 0 M c u n 


•x 1A yW t_ 


J + x + ^ e air + (1 + e air) 




(27) 


ri*CO = PCO^exh ( 28 ) 

Pno -/(PLR - PLRnom) + pR°o m (29) 

fl/vo = PNO Vexh (30) 


4.4. Exhaust gas heat losses model 

For the model, different temperatures are required: the fuel 
temperature, the combustion air temperature, exhaust gas tem¬ 
perature and ash temperature. The fuel and air temperature are 


^h 2 o 


rfifuei 


+ (1 - w) 


y Mh 2 o \ 

2 Mc x \\yO z ) 


( 20 ) 


CO 2 , N 2 and O 2 mass flows in exhaust gas are calculated by using 
following Eqs. (21)—(23): 


™co 2 =^1(1 -w) 


xM, 


CO, 


M 


C X HyO Z 


( 21 ) 


m N 2 


m o 2 


nifuei(l +e a ir)(l -w)(x + | 

nif U eieair(l - w )( x + |-f) 


z\ Xn 2 Mn 2 


2 / xo 2 M C x H y O z 
Mp 2 

M C X HyO Z 


( 22 ) 

(23) 


The combustion air excess which is required to ensure a good 
quality of combustion is often unknown. To reach this parameter, 
Eq. (24) gives the link between this air excess and the oxygen rate in 
the wet exhaust gas. This oxygen rate can be easily reached by using 
an exhaust gas analyser. 


Po 2 


[(*+3-i)£+*+{] 

W M C x H y O z 

(1—w) Mh 2 o 

(*+!-§) 

t 1 P 4 1+ S)] 



(24) 


The oxygen rate in the exhaust gas is controlled by the device by 
adjusting the fan speed. For the tested pCHP device, the oxygen rate 
in wet exhaust gas is equal on average to 6.2% for any PLR that leads 
to an air excess of 48%. The ash and exhaust gas mass flows are 
calculated by Eqs. (25) and (26). For the calculation of the exhaust 
gas mass flow, the ash mass flow is neglected. 

^ash = Xash ^fuel (25) 


^exh ^fuel + tn a j r 


(26) 


Table 4 

Parameters of the model. 



Parameter 

Value 

Unit 

Steady state parameters 



1 

PLRnom 

100 

% 

2 

m nom 

'“fuel 

3.894 

kg/h 

3 

• nom 

Qhx 

16,559 

Wa 

4 

nnom 

1 gross 

1707 

W e l 

5 

pnom 

1 aux 

174 

W e l 

6 

-pnom 

1 exh 

168 

°c 

7 

e air 

0.47 

— 

8 

a 

0.03483 

(kg/h)/100 

9 

b 

152.18 

Wth/100 

10 

c 

16.31 

Wel/100 

11 

d 

0.784 

W e l/100 

12 

e 

0.605 

°C/100 

13 

f 

0.79 

Pg/(n)m 3 /100 

14 

PCO 

36 

t-ig/(n)m 3 

15 

„nom 

^NO 

178 

Pg/(n)m 3 

Transient : 

state parameters 



16 

PLRstart 

27 

0 / 

Zo 

17 

PLRstop 

0 

0 / 

Zo 

18 

PLPmodu 

8 

% 

19 

Qhx 

6 

W t h/s 

20 

P' 

1 gross 

0.6 

Wel/S 

21 

A fQ 
^ L start 

450 

s 

22 

a t-P 
^ l start 

670 

s 

23 

AfQ 

m stop 

2000 

s 

24 

A fP 
^ l stop 

200 

s 

25 

'start 

30 

s 

26 

T P 

'start 

8 

s 

27 

' stop 

800 

s 

28 

k p 

,v stop 

3 

W e ,/s 

29 

E p 

c start 

142 

Wh e l 

30 

E p 

'-"stop 

38 

Wh e l 

31 

mvgj rt 

3.4 

kg/((n)m 3 /s) 

32 

my NO rt 

0.20 

kg/((n)m 3 /s) 

33 

K“ p 

3.7 

kg/((n)m 3 /s) 

34 


0.16 

kg/((n)m 3 /s) 
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defined freely by the user to constant or variable values. The 
exhaust gas temperature is assessed by a linear relation (Eq. (31)) 
depending on PLR. 

^exh - ^exT + e( PLR - PLRnom) (31) 

There are different energy flows created by the incoming and 
out coming components of the system. Here, these energy flows are 
called the enthalpy flows and they are linked to a temperature 
difference between the combustion reactants (fuel and air) and the 
combustion products (ash, CO 2 , N 2 , O 2 and CO). These enthalpy 
flows are calculated by using specific heats or directly specific en¬ 
thalpies. To determine the enthalpy flows of the biomass fuel, 
Wenzl 34] gives an empirical formula of the specific heat of dry 
and wet biomass fuel: 


25,000 

20,000 

^ 15,000 

u 

4) 

£ 

® 10,000 
5,000 

0 



□ Hexhlat 

□ Hdotexh 

■ Qloss 
EQHX 
OPnet 

■ Paux 


PLR [%] 


Fig. 11. Energy balance of the model of the WP pCHP device. 


4.5. Validation of the model and discussion 


Cfuel = (l H4 + 4.86 T fuel )(l - w) +4180w (32) 

The enthalpy flow of the biomass fuel becomes: 

P^fuel = ^fuel Cfuel (^fuel ~ ^ref) (33) 

The enthalpy flow of the latent heat losses is calculated by using 
latent heat of water: 


Hexh = rh^oU (34) 

The exhaust heat losses linked to the others components are 
assessed by using the specific enthalpy of each combustion product. 
These specific enthalpies are extracted from JANAF tables [35] 
which give the specific enthalpy of each component i according 
to its temperature in K. The reference temperature is taken at 273 K. 

" ,= ^( H,(r, ) _H, ( Tref )) (35) 

The exhaust gas heat losses are not recovered here that's why 
this quantity is not included in the thermal production however the 
assessment of this quantity gives the possibility to integrate an 
additional system to recover this heat and to increase the efficiency 
of this system. 

The heat losses by the skin of the pCHP device are worked out by 
Eq. (36). The values have been compared to basic heat losses cal¬ 
culations using a heat exchange coefficient and a temperature dif¬ 
ference between the mean internal air temperature and the 
ambient temperature. By default, this quantity is not recovered by 
considering the device is not located in the heated space, but there 
is the possibility to use this value as internal gain in a building 
dynamic thermal simulation. 


Q-loss — ^fuel T- H a j r + Pf ue i Pgross Q.HX H ex h (36) 

Finally, the model has 34 parameters, 15 for the steady state 
characterization and 19 for the transient states characterization. 
Table 4 sums up all the 34 parameters with the values. 


Table 5 

Maximum deviations between the model and the experimental results in steady 
state. 


PLR [%] 

22 

30 

41 

50 

61 

70 

80 

90 

98 

Efuel [%] 

-0.6 

-1.2 

0.9 

0.1 

-3.2 

-3.2 

-0.6 

1.8 

0.1 

Eth [%] 

-0.2 

-0.6 

-0.3 

2.5 

-3.6 

-2.1 

-0.2 

2.6 

1.0 

Eel m 

3.5 

0.2 

3.0 

-2.7 

-3.1 

-3.6 

-1.2 

3.6 

1.6 

Tex h m 

6.0 

-3.3 

-2.7 

-1.5 

1.0 

-1.4 

1.2 

5.2 

3.0 


The performances of the WP pCHP have been compared in 
steady states to the model. Maximum deviations of 3.2% has been 
found for the fuel energy, 3.6% for the thermal energy, 3.6% for the 
gross electrical energy and 6% for the exhaust gas temperature 
( ^able 5). The minimum duration tests were 3 h. 

The model has been compared to dynamic experimental tests 
which were not used to calibrate the model. The thermal and 
electrical energies produced on the whole cycle are compared to 
the model for 2 cases. Two characteristic test durations have been 
used: the minimum duration test to get good cycle electrical effi¬ 
ciencies (Fig. 12) and the maximum duration test in order to fully 
load the heat storage tank. The first case is a 3 h long test at 70% of 
PLR and the second is a 9 h long test at 70% of PLR. The model 
presents an acceptable accuracy with a maximum deviation of 5.3% 
for the electrical energy of the test 1 (Table 6). 

The R-square value (determination coefficient) of each trend line 
for the steady and transient states is given on Figs. 3 and 4 and 6-9. 

Fig. 11 shows the energy balance of the pCHP device for different 
PLR. The main heat losses are the exhaust gas heat losses. A heat 
recovery on exhaust gas by using a condensation boiler could 
improve the global efficiency of the system. The heat production is 
preponderant compared to the electrical production which is low 
that is mainly linked to the used technology (steam engine) and to 
the power scale. The heat losses are low that let to reach high ef¬ 
ficiencies even at low PLR (Fig. 5). 

The initial motivation is to provide a dynamic model to be 
suitable with annual dynamic thermal simulations for building 
energy simulations and to correctly model the behaviour of a WP 
pCHP to cover heat and power needs with power needs which can 
be very variable. Fig. 12 shows the ratio between the electrical ef¬ 
ficiency of an operating cycle and the nominal electrical efficiency 
according to the cycle duration. This normalized efficiency takes 
into account that the electrical energy which is consumed during 
the start-up phase (starter, HP pump) has to be recovered during 
the beginning of the engine operation before considering that the 
electrical production supplies the building. The cycle should last at 



Fig. 12. Normalized electrical efficiency of a cycle according to the operating time of 
the cycle. 
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Table 6 

Comparison between the model and experimental tests on 2 dynamic cycles. 



£P xp [kWh e ,] 

£ ldell kWh ell 

Deviation [%] 

£g, p [kWh th ] 

CdeilKWhthl 

Deviation [%] 

Test 1 

2.66 

2.80 

5.3 

31.32 

32.60 

4.1 

Test 2 

10.69 

10.65 

-0.4 

109.93 

107.10 

-2.7 


least 3 h to obtain interesting electrical efficiencies (90% of the 
nominal value, Fig. 12). The delay time before the engine triggering 
being at least 600 s, the modulation and a heat storage system 
appear as essential to take maximum advantage of these systems. 

5. Conclusion 

A WP pCHP device has been tested in order to know its ener¬ 
getic and environmental performances. Then a data driven model 
of a steam engine WP pCHP, based on previous experimental 
studies, has been developed and implemented in the TRNSYS 
numerical environment. Both the steady state and the transient 
phases are modelled according to the PLR. It has been shown that 
the inlet water temperature and the mass flow rate of the cooling 
water have an insignificant influence on the energetic and envi¬ 
ronmental performances of the device. Experimental in¬ 
vestigations lead to an electrical efficiency of about 9% and a global 
efficiency of about 95% at full load. The PLR is the main parameter 
having a significant influence on the performances of the device. 
The electrical, thermal and global efficiencies are lower at the 
minimum PLR with an electrical efficiency of about 7.2% and a 
global efficiency of 88%. The pollutant emissions (CO, CO 2 , NO and 
ash) have also been studied. The model just needs 34 experi¬ 
mentally reachable parameters to well characterize the steady 
state and the transient phases by assessing all the energy flows 
involved in the system and the pollutant emissions. The study of 
the transient phases highlights that the cycle durations have to be 
long to get high electrical efficiencies and low pollutant emissions 
impacts. The modelling of the transient phases will let to imple¬ 
ment more realistic numerical studies in order to assess the 
relevance of these systems when they are coupled with buildings. 
This model will be used to carry out optimization studies in order 
to optimize the coupling between buildings and this kind of 
innovative device by considering energetic, economic and envi¬ 
ronment criteria. The main goal is to optimize both the power and 
the heat productions and improve the power self-consumption 
while ensuring the thermal comfort. 
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Nomenclature 

c specific heat (J/kg/K) 

E energy (Wh) 

e a j r air excess (-) 

H enthalpy flux (W) 

H specific enthalpy (J/mol) 

LHV lower heating value (kWh/kgf ue i) 

L v water latent heat (J/kg) 

M molar mass (kg/mol) 

m mass (g) 

mv mass to volume flow rate ratio (kg/((n)m 3 /s) 


ifl 

mass flow (kg/s) 

P 

power (W) 

P 

power ramp rate (W/s) 

p o 2 

oxygen volume rate (%) 

PLR 

part load ratio (%) 

Q 

thermal power (W) 

&' 

thermal power ramp rate (W/s) 

t 

time (s) 

T 

temperature (°C or K) 

V 

volume (m 3 or m 3 /kg) 

V 

volume flow (m 3 /s) 

V m 

molar volume (m 3 /(n)mol) 

w 

water content (kg W ater/kgf U ei) 

X 

molar fraction of carbon (-) 

Y 

mass fraction (-) 

y 

molar fraction of hydrogen (-) 

z 

molar fraction of oxygen (-) 

Gree1<s symbols 

At 

time (s) 

V 

efficiency (%) 

P 

density (kg/m 3 or pg/(n)m 3 ) 

a 

power to heat ratio (-) 

T 

time constant (s) 

X 

volume fraction (-) 

Subscripts and exponents 

air 

air 

amb 

ambiant 

ash 

ash 

aux 

auxiliary 

cw 

cooling water 

dry 

dry part of the biomass fuel 

el 

electrical 

exh 

exhaust gas 

exp 

experimental 

fuel 

fuel 

g 

global 

gross 

gross part of the electrical production 

HX 

heat exchanger 

i 

inlet 

lat 

latent 

loss 

losses 

model 

modelling 

modu 

modulation 

net 

net part of the electrical production 

nom 

nominal 

0 

outlet 

P 

electrical power 

PV 

process value 

a 

thermal power 

ref 

reference 

SP 

set point 

start 

start phase 

stop 

stop phase 

tank 

tank 

th 

thermal 
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Appendix 1. Steady state to transient states transition conditions. 



Current mode 

Future mode 

Trigger 

Stop 

Start-upl 

tPV ^ T SP 
tank ^ tank 

Start-upl 

Start-up2 

End for t - t s tart > At s ^ art 

Start-up2 

Start-up3 

End for t - t star t > At£ art 

Start-up3 

modulation 

End for t - t stan > 5r£ art 

End for Q HX = Q HX 

Modulation 

Steady state 

Steady state 

Cooling 1 

End for t - t stop > Atf top 

Cooling 1 

Cooling 2 

End for t - t stop > At s ^ op 

Cooling 2 

stop 

End for tstart — fstop > At^ op 


Appendix 2. Algorithm chart 

- internal parameters surrounded in dotted orange line. 

- 34 external parameters to be identified surrounded in solid or¬ 
ange line. 

- main outputs surrounded in blue line. 

- main inputs surrounded in black line. 
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inputs 


PLR 


m, 


T 

1 CW,l 


fuel 


T 

1 air 


Y c 

Y h 

Y 0 

Yash 

W 


outputs 





Cooling 

phases 


J PLR„ nm , b , 

'w 





PI R r p nom 

r Lii\ norn , l , i aross 

-* 


>! £CW I" 


gros 
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-> 

■> 
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1 Li£Y nom , M 
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nom> 1 exh 1 ~ 14 
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m CO 
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Yn-exh - >| JANAF 


tables 

LHV -> m fuei > e air c air |->- 

> 





M cw,o 

Qhx 

p 
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p 

1 net 

rn-co 

rn N0 


rjsens 

n exh 


fuel 


is 

n air 

rh C02 

frlat 
n exh 


H 


fuel 


PLRnom> PLPstart> a > PJuel 


PLR nom , pLR start , b, Q n „T, T 


start’ 


pi d pi nstart r pnom T P 
r L,i\ n 0 m> r L,r ' load > r gross’ 1 start > LSl -start 


gross 


PIH PIP rl p nom p" Ar 
' L,n nomi r L,n start’ ‘ aux > E start > LAL start 


mstart a *P 
ln CO i nL start 


mstart a t P 
m N0 , Al start 


* fuel 

Qhx 

p 

1 gross 

Pnet 

rn-co 

ih N0 


PLRnom’ PLRmodu’b, QIT’Q'hX 


PI R r p nom p' 
r Lii\ n 0 m’ c » r grossf r gross 


Qhx 


gross 


PLR nom ,PLR stop , b, Q n H T, 4 op , A t Q stop 


> pjp PI R r p nom t" 

r r L,i\ s top^^r gross’ n -stop 


PI R PI R rl pnom pP \ pe \ t P 
r L,n nom> r L,n stop’ **< r aux > "stop > LXl -stop’ LiL start 


> 

-> 


m p A t p A 

,n CO ’ ^ L starV Al stop 


Stop . p . Q 
rn NO > Alstart> Al stop 



At first, according to the control strategy (heat led, power led), 
the PLR set point is calculated: 


PLR 


SP 


Qhx - QHx m 


+ PLR 
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Appendix 3. Modelling representation of the start-up phase and the cooling phase. 
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